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ABSTRACT 

We present a method of analysing the correlated X-ray and optical/UV variability in 
X-ray binaries, using the observed time delays between the X-ray driving lightcurves 
and their reprocessed optical echoes. This allows us to determine the distribution 
of reprocessing sites within the binary. We model the time-delay transfer functions 
by simulating the distribution of reprocessing regions, using geometrical and binary 
parameters. We construct best-fit time-delay transfer functions, showing the regions 
in the binary responsible for the reprocessing of X-rays. We have applied this model to 
observations of the Soft X-ray Transient, GRO J1655-40. We find the optical variability 
lags the X-ray variability with a mean time delay of 19.3±2.2 seconds. This means 
that the outer regions of the accretion disc are the dominant reprocessing site in this 
system. On fitting the data to a simple geometric model, we derive a best-fit disk half- 
opening angl e of 13.5^9 « degrees, wh ich is similar to that observed after the previous 



outburst by prosz fc Bailyn (1997) . This disk thickening has the effect of almost 
entirely shielding the companion star from irradiation at this stage of the outburst. 

Key words: accretion, accretion discs - binaries: close - stars: individual: Nova Sco 
1994 (GRO J1655-40) - ultraviolet: stars - X-rays: stars 



1 INTRODUCTION 

X-ray binaries (XRBs) are close binaries that contain a rela- 
tively un-evolved donor star and a neutron star or black hole 
that is thought to be accreting material through Roche-lobe 
overflow. Material passing through the inner Lagrangian 
point moves along a ballistic trajectory until impacting onto 
the outer regions of an accretion disk. This material spirals 
through the disk, losing angular momentum, until it accretes 
onto the central compact object. X-rays are emitted from 
inner disk regions via thermal bremsstrahlung with an ef- 
fective temperature ~ 10 8 K. The X-ray flux depends on the 
mass transfer rate, which in turn depends on the structure 
of the disk and it's ability to transport angular momentum. 
In XRBs the structure of the accretion disc is governed by 
irradiation. 

Much of the optical emission in XRBs arises from repro- 
cessing of X-rays by material in regions around the central 



compadt object. The disk is highly ionized and out-shines 
the donor star. Light travel times withm the system are ot 
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order 10s of seconds. Optical variability may thus be de- 
layed in time relative to the X-ray driving variability by 
an amount characteristic of the position of the reprocessing 
region in the binary and the geometry of the binary. The 
optical emission may be modeled as a convolution of the 
lightcurve of the X-ray emission with a time-delay transfer 
function. 

This time delay is the basis of an indirect imaging tech- 
nique, known as echo tomography, to probe the structure 
of accretion flows on scales that cannot be imaged directly, 
even with current interferometric techniques. Echo mapping 
has already been developed to interpret lightcurves of Ac- 
tive Galactic Nuclei (AGN), where time delays are used to 
resolve photoionized emission-line regions near the compact 
variable source of ionizing radiation in the nucleus. In AGN 
the timescale of detectable variations is days to weeks, g iving 



a resolution in the transfer functions of 1-10 light days (Kro- 



lik et al. (1991); Home, Welsh fc Peterson (1991). In XRBs 



the binary separation is ~ light seconds rather than light 
days, requiring high-speed optical and X-ray lightcurves to 
probe the components of the binary in detail. The detectable 
X-ray and optical variations in the lightcurves of such sys- 
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terns are also suitably fast. Extreme examples of this rapid 
va riability arc X-ray bursts, for example from Cygnus X- 
2 (Kuulkers, van der Klis & van Paradijs 1995), where the 
X-ray flux can increase by > 50% with a rise time of 2-3 
seconds and a duration of ~5 seconds. Time-delayed optical 
bur sts have been seen clearly in the object 4U/MXB 1636 



53 (Pedersen et al. (1982); Lawrence et al. (1983) 



Matsuoka 



et al. (1984)). Pulsed X-ray emission from Hercules X-l pro- 



duces faint optical pulsations that ar e thought to be echoes 



from the irradiated companion star ( Middleditch 

i97ey 



Nelson 



Recently, Hynes et al. (1998b) found correlated time- 
delayed X-ray and UV variability in the lightcurves of the 
Soft X- ray transient GRO J1655-40, using RXTE and HST. 



the compact object. These photoionize and heat the sur- 
rounding regions of gas, which later recombine and cool, 
producing lower energy photons. Hard X-rays that penetrate 
below the photosphere emerge as continuum photons, with 
an energy distribution characteristic of the temperature of 
the photosphere. Soft X-rays that are absorbed above the 
photosphere emerge as emission line photons from a tem- 
perature inversion layer near the surface of the disk. Such 
irradiation will have the effect of changing the overall form 
of the orbital lightcurve as the aspect of the hot, irradiated 
regions changes with the binary orbit, this effect is very no- 
ticeable with the changing irradiation of the inner face of 
the companion star in Hercules X-l, during the 35-day pre 



cessio n phase of the tilted accretion disc (Boynton et al 



The da ;a, centred around binary phase 0.4, shows a mean 
time-delay of 14 seconds and an RMS delay of 10.5 seconds. 
The mean time delay is consistent with reprocessing in the 
outer regions of the accretion disk, while the relatively large 
range of delays suggests that there is relatively little com- 
panion star reprocessing. This may occur if a thick outer 
accretion disk shields the companion star from the X-ray 
source. 

In this paper we present a simple geometric model for 
the time-delay transfer functions of XRBs, using a synthetic 
binary code. We analyse correlated X-ray and UV variability 
in GRO J1655-40, using our computed transfer functions, 
to constrain the size, thickness and geometric shape of the 
accretion disk. 

Section || describes the features seen in typical optical 
lightcurves of LMXBs and how these affect our analysis. Sec- 
tion ^ describes the reprocessing of X-rays in XRBs in more 
detail, including the reasons for the time delays between the 
X-ray and optical variability. Section ^| describes the model 
created to describe the time delays found, while Section ^ 
describes an alternative method that uses Gaussian transfer 



1972 ) . While these features are all vital in interpreting the 



long timescale orbital lightcurves of X-ray binaries, many 
of them are not important in this work, which deals with a 
very small range of binary phases (less than 1%). 

The short duration of the observations is also important 
when considering the effects of other transient features, such 
as star spots, whose lifetime is much longer than the dura- 
tion of our observation and whilst it will undoubtedly make 
a small effect on our model fits, the uncertainties involved in 
including them far outweigh any benefits to the model from 
including them. 

In our model, we have assumed that all the optical 
variability is in fact caused by the reprocessing of X-ray 
irradiation. Furthermore, we have assumed that the ther- 
mal component will dominate the reprocessed emission in 
our broadband observations. These assumptions are clearly 
not strictly true, observations of cataclysmic variables and 
other interacting binaries, where irradiation is no longer 
the dominant source of optical emission, also show non- 
periodic optical variability superimposed on the periodic 



orbital lightcurves (e g. the Nova-like AE Aquarii, Welsh 



is followed by a discussion of the results in Section |^. 



2 OPTICAL LIGHTCURVES 

The lightcurves of X-ray binaries contain many temporal 



functio: is. Section |6 shows the results of this analysis. This Horne fc Oke (1993) ). However, while such non-correlated 



and spectral features (see van Paradijs & McClintock (1995) 



for a review of many of these), some periodi c and others 
purely random or quasi-periodic in nature (See van der Klis 



(2000) for a review). The periodic features, such as eclipses 



and the minima and maxima of lightcurves are due to the 
orbital motion of the components of the binary around the 
centre of mass of the system. These periodic features are 
relatively easy to model and much success has been made 
of such models in determining the geometric parameters of 
binaries using the information in multi-colour lightcurves. 
Refinements have been made to include effects such as limb 
darkening, caused by an decreasing source function with pro- 
jected radius from the stellar core, and gravity darkening, 
which causes a change in temper ature of the distorted sh ape 
of the roche-lobe filling star (eg. Orosz & Bailyn (1997)). 



In the case of X-ray binaries the intense X-ray flux from 
the regions surrounding the compact object also affect both 
the periodic and non-periodic features of the observed op- 
tical lightcurves. In the standard model of reprocessing, X- 
rays are emitted by material in the deep potential well of 



features will introduce noise into out analysis, thus increas- 
ing the absolute value of the badness-of-fit of our models, 
they will have little effect on the relative values, which have 
contributions from all points in the lightcurve. Similarly, the 
optically thin emission is only a small fraction of the total 
broadband emission and for this reason is assumed to vary 
in phase with the thermal emission. 



3 REPROCESSING OF X-RAYS 

The reprocessed, optical emission seen by a distant observer 
is delayed in time of arrival relative to the X-rays by two 
mechanisms. The first is a finite reprocessing time for the X- 
ray photons and the second is the light travel times between 
the X-ray source and the reprocessing sites within the binary 
system. 



3.1 reprocessing times 

The average reprocessing time fo r line photons is given b y 
the average recombination time, (Hummer & Seaton 1963) 
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In the accretion disk the high electron density, n e ~ 
10 15 cm -3 , ensures that the time-scale for reprocessing of 
line photons is short compared to the overall time delay. 

The continuum photons which scatter from deeper 
within the accretion disk undergo a 'random walk' before 
escaping through the photosphere, leading to a longer re- 
processing time. The exact determination of this continuum 
reprocessing time is complicated, requiring detailed model 
atmosphere calculations, which is outside the scope of this 
paper. Howeve r, the short time de lays between X-ray and 
optical bursts ( Pedersen et al. 1982 ) and the observ ations of 
1.24 sec ond optical pulsations from Hercules X-l (Chester 
197S |), | mply that a significant fraction of the reprocessed 
optical photons emerge from the reprocessing site within 
~ 0.6 seconds of the absorption of the incident X-ray pho- 
tons. This delay is smaller than the measured uncertainty 
in the mean delay for the systems and so we have treated 
the reprocessing as instantaneous. We have also treated the 
reprocessing of X-rays as "passive" reprocessing, where the 
absorbed X-rays do not affect the structure of the material 
in the binary. 



3.2 light travel times 

The light travel times arise from the time of flight differences 
for photons that are observed directly and those that are re- 
processed and re-emitted before traveling to the observer. 
These delays can be up to twice the binary separation, ob- 
tained from Kepler's third law, 



- = 9.76s 
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where a is the binary separation, M x and M c are the masses 
of the compact object and donor star, P is the orbital period. 
In LMXBs the binary separation is of the order of several 
light seconds. 

The time delay r at binary phase <j> for a reprocessing 
site with cylindrical coordinates (R, 6, Z) is 



t(x, <p) — (1 + sm i cos( 



6)) cos i 



(3) 



where i is the inclination of the system and c is the speed of 
light. This can also be expressed using the position vector, 
x and the unit vector, e((J>), pointing toward the earth, 

e((/>).x 



"(x,0) = 



(4) 



The X-ray driving lightcurve, fx(t), is described as the 
sum of a constant and a variable component, 



/»(t) = /- + A/«(t). 



(•») 



The reprocessed lightcurve, f v (t), is similarly divided into 
two components. The relationship between the two is given 
by 



(6) 



where ^ v (A, r, <j)) is the time delay transfer function. 
This transfer function is the strength of the reprocessed vari- 
ability delayed by r relative to the X-ray variability. 

The dynamic response function is found by considering 
how a change in X-ray flux drives a change in the reprocessed 



flux. We can define the dynamic time delay transfer function 
to be, 



W„(A,x,A/«(t-T)) 



dfl(x, cf>) . S(t - r(x, 0)) 



(7) 



where r(x, 4>) is the geometric time delay of a reprocessing 
site at position x, see Equation [| In the next section we 
describe the model X-ray binary code we have been using 
and how we have used this to find the reprocessed flux. 



4 MODEL X-RAY BINARY CODE 

We have developed a code to model time delay transfer func- 
tions based on determining the contributions from different 
regions in the binary. In this section we describe the mod- 
els used to construct the individual regions of the binary; 
the donor star, the accretion stream and the accretion disk. 
The code uses distances scaled to the binary separation in 
a right-handed Cartesian coordinate system corotating with 
the binary. The X-direction is along the line of centres for 
the binary, the Y-direction is perpendicular to this in the 
orbital plane of the binary, so that the X-ray source is at 
(0,0,0) and the centre of mass of the donor star is at (1,0,0). 
Each surface element is a triangles, characterized by its area 
dA, orientation n, position x and temperature T. 

We calculate the total monochromatic flux by summing 
up contributions from all visible elements. The area dA and 
the normal vector n for the triangular panel are calculated 
and then the projected area of the panel is calculated using 
the projected earth vector e(<p,i). The effects of occultations 
by regions in the binary are also considered. Therefore the 
observability, 0(x, cj>), is given by, 



0(x,<p) = dA(x).e(<£) 



(8) 



which is the foreshortened area of pixel x, observable at 
phase <j>. This is related to the solid angle of the observed 
pixel, dfl by the relation 



df2(x, 



0(x, 



D 2 



(9) 



where D is the distance to the source. If dQ,(x, (f>) < then 
the panel is not visible to the observer and does not con- 
tribute to the flux. The monochromatic intensity is calcu- 
lated using the Planck function, 

B ^ T ) = T7T \ 71 ( 10 ) 



^ 3 [«p(j&)-l] 



This is scaled using the projected area of the panel |0(x, <f))\, 
as seen by the observer. The standard linear limb-darkening 
law is assumed, 



L(u, a) 



- u + u cos (a) 
l + u/3 



(11) 



where u is the linear limb-darkening coefficient, assumed to 
be 0.6, and a is the angle between the normal vector n and 
the earth vector §(</>, i), 

cos a — n.e. (12) 

The response curve for a given detector, P (A), is com- 
bined with the limb-darkened Planck function to create the 
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synthetic reprocessed flux from each visible triangular el- 
ement. The total detected flux due to reprocessed X-rays 
from a single panel, F v (X,t), is given by, 



F U (X,T) = / B v (X,T) P (X) L (u, a) dQ (x, 0) dX, 



(13) 



where L (u, a) and dQ, (x, </>) are the expressions for the 
limb-darkenin g a nd the solid angle of the exposed panel, 
see equations [ll] and ^ respectively. In the next sections 
we describe the geometric model used to calculate the time 
delay for a given panel and the effects of irradiation which 
are used to determine the contribution of each panel to the 
final transfer function. 



4.1 Donor Star 

The donor star is modeled assuming it fills its critical Roche 
potential, so that mass transfer occurs via Roche lobe over- 
flow through the inner Lagrangian point. Optically thick 
panels are placed over the surface of the Roche potential. 
The panels are triangular so that the curved surfaces of the 
binary are map ped more accurately tha n is possible using 
4-sided shapes (Rutten & Dhillon 1994). These panels are 



equally spaced in longitude and latitude across the surface 
of the star. In order to correctly specify the temperature of 
each panel on the face of the tidally distorted star, one must 
take into account the degree o f gravity darkening. Using von 
Zeipel's theorem (Zeipel 1924), for the relationship between 
the local gravity and the local emergent flux one finds that 
the relationship between the local temperature, T e , and the 
gravity, g, is 



T e 4 (x) « g(x) 



(14) 



As a consequence the temperature at any point on the star 
is given by 

(15) 



T(x) 



Tpole 
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where T po ; e and g po ie are the temperature and gravity of 
the pole of the star. The "gravity darkening exponent" 



is 0.25 for stars with fully radiative envelopes (Zeipel 1924 



as is the case in our models) an d 0.08 for stars with fully 
convective envelopes (Lucy 1966). T po i e is taken to be the 
effective temperature of a field star with a similar spectral 
type to the donor star. 



4.2 Accretion stream 

The accretion stream is modeled by following the ballis- 
tic trajectories of 4 test particles. The thickness (w) of 
the stream defines the initial positions of the test particles. 
These test particles determine the 'width' of the stream (Its 
deviation from the line of centres of the binary, in the plane 
of the binary, or y-direction) and the 'height' of the stream 
(it's extent in the direction normal to the plane of the binary, 
the z-direction), assuming the stream is symmetric about the 
x-y plane. 

The particles start at the LI point with a small ve- 
locity in the direction of the compact object (-v,0,0), 
from positions (R(L1),0,0), (R(Ll),w,0), (R(Ll),-w,0) and 
(R(Ll),0,w). The trajectory is cut into discrete steps, with 



the step size as a parameter of the code. The velocity and 
position of each particle are determined from the Roche po- 
tential after each step. The stream is curtailed when one of 
two criteria are reached; (1) the stream has collapsed verti- 
cally, or (2) the core trajectory is moving outwards, ie. the 
stream has passed the compact object without collapsing 
vertically. 

The unirradiated accretion stream is assumed to have 
a constant temperature T s along it's length and the effects 
of irradiation are considered in the same way as those of the 
donor star in Section i.4. 



4.3 Accretion disk 

The disk thickness is assumed to increase with radius from 
at R = Ki n to H «i at R = Rout, with the form, 



H = Rout (H\ / R-R in \ 

V R ) out V Rout — Rin ) 



(16) 



where the parameters are the inner and outer disk radii, 
Ri„ and R ou t in units of R(L1), the half thickness of the 
outer disk (H/R) ut and the exponent which describes 
the overall shape of the disk. The temperature structure of 
the un-irradiated disk is that of a steady state disk, in the 
absence of irradiation, 



Tdisk(R) — Tout 



R 

Rou 



(17) 



where T ou t is the temperature at the outer disk and T; n is 
the temperature of the inner disk. 

The disk is divided radially and azimuthally, into Nr 
and No sections. The monochromatic intensity is again 
calculated using the Planck function, corrected for limb- 
darkening using a linear limb-darkening law, with a constant 
coefficient. The intensity is again scaled using the projected 
area of the panel, for the given values of binary phase (j> and 
inclination i. 



4.4 Irradiation model 

The effective temperature of a region at a distance R from 
the X-ray source, assumed in our model to be a point source 
located at the centre of the accretion disk, is found from the 
accretion luminosity for a typical LMXB, 



and 



Ml- A) 
4naR 2 



GM X M 



Rn 



(18) 



(19) 



where T x is the temperature, A is the albedo, r\ is the 
efficiency, M x , the mass of the compact object, M the ac- 
cretion rate onto the compact object, R ns is the size of the 
compact object and R is the distance between the compact 
object and the irradiated element. This is normalised using 
the binary separation, a, the distance between the centres of 
mass of the stars, as is the coordinate system for the binary. 
In the case of Scorpius X-l, this gives T x ~ 10 K for a 1.4 
Mq neutron star (Rns ~ lOfcra) accreting 10 _9 MQt/r _1 , 
with an efficiency r\ — 0.1, an albedo of 0.5, at a distance 
equal to the binary separation of 3.4 x 10 11 cm. 
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The irradiation of the binary takes place in three stages. 
The first stage is to calculate the temperature structure of 
the binary in the absence of any irradiation. This is done 
with characteristic temperatures for the donor star (from 
it's spectral type) and the accretion stream and disk. The 
temperature structure of the disk is assumed to that for an 
unirradiated disk as given in equation [l?]. The surface ele- 
ments of the binary exposed to X-rays are determined by 
projecting the binary surfaces onto a spherical polar repre- 
sentation of the sky, as it appears from the X-ray source. 
Each triangular element is mapped to the sky starting with 
the one furthest from the source and ending with the trian- 
gle closest. Those elements remaining visible and unocculted 
on the sky map are irradiated. The change in effective tem- 
perature of an element is scaled by the projected area with 
respect to the X-ray source at a distance R from the source. 
Hence the temperature after irradiation is given by, 



T 4 = T*cos0f-! 

it 



<-'// 



(20) 



where T is the temperature of the panel, 6 X the angle be- 
tween the line of sight from the central source and the nor- 
mal to the surface of the element and T e f / is the unirradiated 
effective temperature of the panel. 

Since we are interested in finding the correlation be- 
tween the variable component of the X-ray and reprocessed 
fluxes, we split the X-ray flux into constant and time- 
dependent components. These components of the flux are 
converted into components of temperature, 



T x (t) = T x + AT x (t) 
where 

A/„(t) 4AT a (t) 



(21) 



(22) 



The second stage is to irradiate the binary with the 
constant component of the X-ray flux. This component of 
the X-ray flux is equated to the mean effective temperature 
of the X-ray source, as given in equation |5o| where T x = 
T x . The third and final stage is to repeat stage two with 
T x = T(t) x , which represents irradiating the binary with a 
time varying component. The difference between stages two 
and three represents the temperature change of the elements 
due to the time varying component of the X-ray flux alone, 
A/„(t). 

Thus the response of a panel to the variable component 
of the irradiating X-ray flux is given by, 

I„(\,x,Af x (t)) = J [B u (X,T x (t)) - B v (\,Tx)] . 

P (A) I (u, a) dQ (x, 0) dX. (23) 

This response is substituted into the expression for the dy- 
namic response given in equation ^. 

4.5 Transfer functions 

In order to transform the reprocessed flux into a time de- 
lay transfer function, we define iso-delay surfaces. These 
surfaces are nested paraboloids around the line of sight to 
the X-ray source, defined by the earth vector, e(<f>,i). The 
parabolic surfaces have a mean time delay r and a width 
St. The mean time delay r in seconds between the directly 




time delay (seconds) 

Figure 2. A plot of time-delay transfer functions as a function 
of binary phase, based on the typical binary parameters for an 
LMXB. The accretion disk has constant time delays in the region 
0-4 seconds, whereas the time delays from the companion star 
are seen to vary sinusoidally with binary phase between 2 and 10 
seconds. 



observed X-ray flux from the central source and the re- 
processed signal from a point with cylindrical coordinates 
(R,0,Z) is 



VR 2 + . . Z . 
T = (1 + sm % cos(<p — 0)) COS I 



(24) 



where i is the inclination of the system, (f> is the binary phase 
and 6 is the angle from the line of centres of the binary. 

We consider models with reprocessing taking place in 
the accretion disk, accretion stream and companion star. 
Each of these regions makes a contribution to the transfer 
function for the system, see Figure |l| for a diagram of the 
results of our code and the calculated transfer functions. 

The phase dependence of r allows us to create time- 
delay transfer functions as a function of binary phase, al- 
lowing us to produce phase-delay diagrams for the system, 
see figure |^. X-rays reprocessed at the companion star have 
a time delay that varies sinusoidally in phase with semi- 
amplitude (a/c) sin i around a mean value a/c. X-rays re- 
processed by a circular disk appear as a phase independent 
contribution to the delay distribution 3>(r, (f>) between the in- 
ner and outer radii of the disk. The accretion stream shows 
up as a non-symmetric contribution that varies roughly si- 
nusoidally with the orbital motion of the companion star 
and can be seen faintly in Figure ^| as the contribution near 
superior conjunction between the outer rim of the disk and 
the inner-face of the companion star. 

The relative intensities of the contributions, represented 
by the area under the transfer function, constrain the geo- 
metric parameters of the system, especially the contribution 
from the accretion disk, which is probably the most impor- 
tant region for reprocessing of X-rays in interacting binaries. 

Thus the shape of the iso-delay surfaces when projected 
onto the plane of the binary depends on the inclination of the 
system. For a face-on disk (i = 0°) the projected iso-delay 
surfaces are simple circles on the disk. As the inclination of 
the system increases these surfaces become elongated along 
the line of sight to the X-ray source, until finally they form 
parabolae if viewed from edge-on (i = 90°). This inclination 
dependence is shown in the transfer functions in Figure pi 
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Figure 1. Left, model X-ray binaries, based on typical binary parameters for an LMXB, showing iso-delay surfaces projected onto the 
irradiated surfaces of the binary. Right, the associated time delay transfer functions, showing the relative contributions from the regions 
highlighted in the model X-ray binaries. 
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Figure 3. Model time-delay transfer functions for an accretion 
disk as a function of binary inclination. The solid line is for a 
binary inclination of 0°, the dashed for 35° and the dot-dash for 
70° 



Delay (sees) 

Figure 4. Model time-delay transfer functions as a function of 
/3, the disk exponent for f}= 1.01(solid line), 1.5(dashed line), 
2(dot-dash line), 2.5(dotted line) and 3(3 dots-dash line) 



The peak in these transfer functions is also a function of 
inclination, with the peak occurring with the time delay of 
the largest iso-delay surface that is contained entirely within 
the disk, 



R 



(1 — sin (i — 5)), 



(25) 



where 5 is the opening angle at the edge of the disk, 8 — 
tan -1 

The disk shape, characterized by the exponent /3, also 
affects the position and shape of the peak in the transfer 
function. The contribution is scaled with the projected area 
of the surface element, therefore as (3 increases the peak 
should move to longer time delay, as the edge of the disk 
becomes steeper and the projected area decreases. This can 
be seen in figure ^| where time-delay transfer functions are 
plotted as a function of (3 with values between 1.01 and 
3. (Note: When /3=1, the surface of the disk is flat, whose 
thickness goes to zero at the position of the compact object, 
therefore the projected area of the disk as seen from the 
compact object is zero, except on the inner edge. For this 
reason f3 is given a value just greater than 1 to show the 
limit as (5 tends towards the case of an un-flared disk.) 



Parameter 


symbol 


unit 


GRO 
J1655-40 


Distance 


D 


kpc 


3.2 


Period 


P 


days 


2.6 


Inclination 


i 


deg. 


69.5 


Mass ratio 


q 




0.3344 


Primary mass 


M :r 




7.02 


Binary sep. 


a 


10 11 cm 


12 


Star Temp 




K 


6500 


Stream Temp 


T s 


K 


5000 


Outer Disk 


Tout 


K 


6000 


Inner Disk 


T 


K 


100000 


Irrad Temp 


T x 


K 


10 5 


Inner radius 


Rin 


R(il) 


10" 4 


Outer radius 


Rout 


R(ii) 


free 


Disk thickness 


H/R 




free 


disk exponent 


& 




free 



Table 1. Parameters used in our model of X-ray Binary, the 



values are taken from 



Orosz fe Bailyn (1997)| 



5 GAUSSIAN TRANSFER FUNCTIONS 

In paper I Gaussian time delay transfer functions were 
used to create synthetic reprocessed lightcurves for the SXT 
GRO J 1655-40. This is a relatively simple form of transfer 
function as it assumes nothing about the geometry of the 
system. These transfer functions have the form, 



*(t) 



2ttAt 



exp 



1 ( T ~ T ° V 

2 V Ar J 



(26) 



Thus the transfer functions have three parameters, the 
mean time delay To, it's variance Ar and the strength of the 
response Vf. These transfer functions were convolved with 
the X-ray driving lightcurve in the same way as the model 
transfer function and the badness of fit between the syn- 
thetic and real reprocessed lightcurves calculated. 

The equivalent parameters ro and Ar can be deter- 
mined for any transfer function, \&(t) by calculating the first 
and second moments. 



to 



Ar 



(r-r ) 2 dr 
f°° 9M dr 

■J — oo v ' 



1/2 



(27) 



(28) 



6 OBSERVATION AND DATA REDUCTION 

The Soft X-ray Transient (SXT) GRO J1655-40 was discov- 
ered in 1994 July when the Burst and Transient Source Ex- 
periment (BATSE) on GRO observed it in outb urst at a 
level of 1. 1 Crab in the 20-200 keV energy band ( Harmon 
et al. 1995). After a period of apparent quiescence from late 



1995 to early 1996, GRO J1655-40 went int o outburst again 
in late 1996 April ( Remillard et al. 1996) ) , and remained 
active until 1997 August. During the early stages of this 
outburst we carried out a series of simultaneous HST and 
RXTE visits. One of the primary goals of this project was 
to search for correlated variability in the two wavebands. 
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The long-term evolution of this outburst argued against sig- 
nificant reprocessing, as th e seemingly anticor related optical 
and X-ray fluxes observed (Hynes et al. 1998a) are not to be 
expected if the optical flux is reprocessed X-rays. Nonethe- 
less, significant s hort term correlatio ns were detected. In our 
previous paper (Hynes et al. 1998b) we have analyzed these 



correlations using both acausal and causal Gaussian trans- 
fer functions. We were able to put constraints on the possi- 
ble regions responsible for X-ray reprocessing. In this paper 
we summarize this previous analysis and develop it further, 
using our synthetic binary model to create more physical 
transfer functions. 



6.1 observations 

The HST and RXTE observations used in this paper were 



first pr esented, along with similar exposures, in ( Hynes et al 



1998b). In paper I we fitted causal and acausal Gaussian 



transfer functions to 4 similar exposures in a similar analysis 
as used in the previous section. The time-delay distribution 
based on all the observations is 14.6 ± 1.4 seconds, with a 
dispersion of 10.5 ± 1.9 seconds. 




Figure 5. The best-fit results for GRO J1655-40 using acausal 
Gaussian transfer functions. Top panel, the normalised X-ray 
driving lightcurve from RXTE. Middle panel, UV lightcurve from 
HST with synthetic UV lightcurve superimposed (thick line). Bot- 
tom panel, residuals of the fit to the UV lightcurve. 



6.1.1 HST 

The HST observations are shown in figure 1(b) of paper 
I, along with a detailed description of the data reduction 
techniques used. In this paper we have used one of the ob- 
s ervations, referred to as exposure 6, shown in Figure 1(b) in 
( Hynes et al. 1998lj ) and reproduced in this paper as figure || 
This observation took place during the June 8 1996 visits, 
using The Faint Object Spectrograph in RAPID mode with 
the PRISM and blue detectors (PRISM/BL), covering the 
spectral range ~ 2000-9000 A. The resulting lightcurve has 
a time resolution of ~ 3 seconds, the absolute time accuracy 
of this data is limited to 0.255 seconds. 



6.1.2 RXTE 

The X-ray data was taken with the PCA onboard RXTE 
on June 8 1996, simultaneous with the HST data. The 
lightcurve was created from the standard-1 EDS mode data, 
using the saextrct task in the FTOOLS software package. This 
mode has a maximum time-resolution of 0.125s, but with no 
spectral resolution. It is desirable for echo-mapping for the 
driving X-ray lightcurve to have a higher time resolution 
than the reprocessed one, so that the time delay transfer 
function can have a time resolution greater than the time 
resolution of the reprocessed lightcurve. A lightcurve with 
a time resolution of Is and an absolute timing accuracy of 
about 8 /xs was extracted. For a full description of the data, 
see paper I. 

6.2 results of modeling 

6.2.1 Gaussian transfer functions 

The results for the fits to the acausal Gaussian transfer func- 
tions are described in depth in paper I and summarized in 
this paper in Table |^, together with their 1-parameter 1- 
sigma confidence regions, in order to compare them to the 
results from our modeling of the X-ray binary. The mean 



Dataset 



Exp. 6 



0.42 
782 



X 2 /(N-3) 


1.230 


TO 


19.3 ± 2.2 


At 




4V1CT 3 





Table 2. Best fit values for the acausal Gaussian time delay 
transfer function fitting to the X-ray and optical lightcurves for 
GRO J1655-40. The label exp. 6 refers to exposure 6 in the origi- 
nal paper, Paper I. 



and rms values for the delay in the Gaussian transfer func- 
tion represent the first and second moments expected for 
the model transfer functions. The best fit Gaussian transfer 
function for expo sure 6 has a mean de lay of 19.3s with an 



rms delay 10.8s. Hynes et al. (1998b) interpreted this re- 
sult as evidence for reprocessing in the outer regions of a 
thick accretion disk that could be thick enough to shield the 
companion star from significant irradiation. If this were not 
the case and we had considerable disk and companion star 
reprocessing we would expect the variance of the transfer 
function to be much higher. 

6.2.2 synthetic binary transfer functions 

The transfer functions created by our X-ray binary model, 



using t he binary parameters determined by Orosz & Bailyn 
(1997) and summarized in Table [l] were used to create syn- 



thetic HST lightcurves by convolving them with the RXTE 
lightcurves. The badness of fit between this synthetic repro- 
cessed lightcurve and the real reprocessed lightcurve, from 
the HST observations was calculated. The fitting again con- 
tained 3 parameters, which were optimized as before. The 
best-fit parameters, together with their 1-parameter 1-sigma 
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Dataset 



Exp. 6 



N 



0.42 
782 



X 2 /(N-4) 1.230 

J^out u - u,j _q 06 

H/R 0.24±gig* 
P 4.7lf 



Table 3. Best fit values for the synthetic X-ray Binary model time 
delay transfer function fitting to the RXTE and HST lightcurves 
for GRO J1655-40. 




Figure 6. A comparison of the best fit transfer functions for 
GRO J1655-40 from our two modeling methods. On the left is 
the synthetic X-ray binary model transfer function and on the 
right is the acausal Gaussian transfer function. 



confidence regions are given in Table |j| We find that the disc 
extends to 67 percent of the distance to the inner Lagrangian 
point, with an opening angle of 14 degrees and an exponent 
of 3.8. 



6. 2. 3 comparison of the two models 

This best-fit solution from our modeling is shown along with 
the best-fit solution from the acausal Gaussian fitting in Fig- 
ure 13. The best- fit solution for this modeling shows a similar 
X 2 to the Gaussian fitting, which shows that our modeling 
of the binary system and simple Gaussian fits are both good 
fits to the data. The first and second moments from the 
model transfer function are, as expected, similar to the best 
fit parameters from the Gaussian fitting. The opening an- 
gle for the disc from our synthetic transfer functions, whilst 
being surprisingly large is comparable with that derived by 
Orosz & Bailyn (1997) from lightcurve fitting shortly after 
the previous outburst in March 1995. They found that the 
opening angle of the disc was 11° (see Table 6 from Orosz & 
Bailyn (1997)) compared to our value of 14°. 



parameters for the binary system GRO J1655-40. These pa- 
rameters are principally the size and shape of the accretion 
disk. This is inferred from the relative contributions to the 
time delay transfer function of the different regions of the bi- 
nary. We have used time delay transfer functions, along with 
the known binary parameters to find best fit solutions to the 
data along with their corresponding confidence regions. 

There is evidence for a larger fraction of disk repro- 
cessing. The geometric parameters determined from our fit 
give an opening angle of 19 degrees and a very flared ge- 
ometric shape, implying that most of the reprocessing is 
taking place in the outer regions of the disc. This means 
that the companion star is almost entirely shielded from X- 
rays, which in turn reduces the mass accretion rate, which 
is driven by irradiation. The observations of GRO J1655-40 
took place during outburst, which could explain the flared 
shape of the outer disk. Another possible interpretation is 
that there is a localized region of enhanced reprocessing in 
the outer disk that is non-axisymmetric that is adding a 
large component from the outer disc region to the trans- 
fer function but is difficult to distinguish from reprocessing 
from a thick disk. This would explain the high value of /3 ob- 
served, as the model attempts to move all the reprocessing 
to the outer disk. The most likely site for this reprocessing 
would be at the disk-stream interaction point, where the rim 
of the disk swells greatly. This inhomogeneity is observed in 
GRO J 1655-40 as X-ray dips around binary phase 0.8. This 
is the site proposed for the reprocessing of the X -ray bursts 
seen in 4U/MXB 1636-53 ( |Matsuoka et al. 1984| ). 

The effect of reprocessing timescales in different regions 
of the binary also needs to be studied in more detail, with 
detailed radiative transfer models. The incidence angle of 
the X-rays to the atmosphere may cause large variations in 
the timescale. Normally incident X-ray photons may be re- 
processed deeper within the companion star and hence take 
longer to diffuse to the surface, than those with a grazing 
incidence angle. This would also affect the ratio of disk to 
companion star reprocessing, as the incidence angles for the 
disk will be predominantly grazing. 

In order to distinguish between these scenarios and de- 
termine the importance of reprocessing timescales it is nec- 
essary to have data with better phase coverage to observe 
a change in the time delay of this enhanced region as a 
function of binary phase and any reprocessing timescale ef- 
fects. It is clear from our analysis that the companion star 
is responsible for a small fraction of the instantaneously re- 
processed flux, see figure [| 

This is the first data showing correlated X-ray and op- 
tical variability with sufficient time resolution and a long 
enough base-line to do this form of echo-tomography. We 
have used the time delayed optical variability observed from 
GRO J1655-40 to constrain the binary parameters of these 
systems and find that both the accretion disk and compan- 
ion star can be important regions for reprocessing in XRBs. 
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